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Abstract 

We study the cosmology of the Brans-Dicke(BD) theory coupled to perfect 
fluid type matter. In our previous works, the case where matter is coming 
from the Ramond-Ramond sector of the string theory was studied. Here we 
study the case where matter is coming from NS-NS sector. Exact solutions 
are found and the cosmology is classified according to the values of 7, the 
parameter of the equation of state and uj, the BD parameter. We find that, 
in string frame, there are solutions without singularity for some ranges of 7 
and UJ. In Einstein frame, however, all solutions are singular. 
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I. INTRODUCTION 



The string theory is beheved to be the most promising candidate to quantum gravity. So 
it is natural to expect that it will resolve some problems inherent in the general relativity 
like the initial singularity problem. In fact, in the regime of Planck length curvature, quan- 
tum fluctuation is very large so that string coupling becomes large and consequently the 
fundamental string degrees of freedom are not weakly coupled good ones. Instead, solitonic 
degrees of freedom like solitonic p-branes [T|] or D-brane [Q] are more important. Therefore 
it is an interesting question to ask whether including these degrees of freedom resolve the 
initial singularity. 

The new gravity theory that can deal with such new degree of freedom should be a 
deformation of standard general relativity so that in a certain limit it should be reduced to 
the standard Einstein theory. The Brans-Dicke theory is a generic deformation of the 
general relativity allowing variable gravity coupling. Therefore whatever is the motivation 
to modify the Einstein theory, the Brans-Dicke theory is the first one to be considered. As 
an example, low energy limit of the string theory contains the Brans-Dicke theory with a 
fine tuned deformation parameter {uj = —1) and it is extensively studied under the name 
of the string cosmology . Without knowing the exact theory of the p-brane cosmology, 
the best guess is that it should be a Brans-Dicke theory with matters. In fact there is some 
evidence for this where it is found that the natural metric that couples to the p-brane 
is the Einstein metric multiplied by certain power of dilaton field. In terms of this new 
metric, the action that gives the p-brane solution becomes Brans-Dicke action with definite 
deformation parameter u depending on p. 

In our previous works ^j^, we studied the gas of solitonic p-brane 0] treated as a 
perfect fluid type matter in a Brans-Dicke theory allowing the equation of state parameter 7 
arbitrary. We had studied the case where the perfect fluid does not couple to the dilaton like 
the matter in the Ramond-Ramond sector of the string theory. Here we study the opposite 
case where matter couples to the dilaton like those coming from NS-NS sector(see reference 
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0] for similar study). Exact solutions are found both in string and Einstein frame and the 
cosmology is classified according to the values of 7 and uj. In string frame we will find 
non-singular solutions for some ranges of 7 and uj. In Einstein frame, however, we will find 
that all solutions are singular, unlike the string frame. 

The rest of this paper is organized as follows. In section II, we construct the action for 
the case when the matter coupled to the dilaton and find analytic solutions for the equations 
motion. In section III, the relation between cosmic time t and parameter r is studied. In 
section IV, we consider the behavior of the scale factor a as a function of r. Using these 
results, in section V, we study the scale factor as a function of the cosmic time t. In section 
VI, we study the asymptotic behavior of a{t) and classify them according to acceleration 
and deceleration phases. Up to the section VI, all analyses were done in string frame. In 
section VII, we investigate the cosmology in Einstein frame. In section VIII, we summarize 
and conclude with some discussion. 



II. ACTION WITH SOLITONIC NS-NS MATTER AND ITS SOLUTIONS 

We begin with the four dimensional Brans-Dicke like string action of which matter is 
coupled to the dilaton field. 

S = [ (fx^^e-^ \r - uj{W(t)f + Lrr\ (1) 



Notice that the action differs from that in reference by the coupling of the matter 
Lagrangian with the dilaton factor. By varying the action, we get equations of motion: 

i?- 2cjVV + ^(V0)^ + L„ = 0. (2) 
Let's choose the metric as 

ds"^ = -Ndt^ + e^""^^^ dxidx' = 1, 2, 3), 



where is a lapse function. We consider perfect fluid type matter whose energy-momentum 
tensor is 

Tn.u = pgfiu + (p + p)UnUv, 

satisfying the conservation law 

P + 3(p + p)q; = 0. 
Using the equation of state, p = jp, we get 

p = poe-^(^+T)". 

So, we can rewrite the action as 
Now we define new variable r by 
Then, the action becomes 



where 



(3) 



= ldr[r,Y" + r,X"-poe-'''l (5) 



^ -6 + 9(1 - 7) + ^(1 - jfcu = ^(1 - j)\u; - curj, 

r2 = 6 + 3(1 - 7)cu = 3(1 -i){u;- a;rj, 
. _^ rl _ 3(2a; + 3) 



-2X = 3(l-7)a-20, 

= (7) 
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From this action we get equations of motion: 



Y" = 0, 
X" - = 0. 



The constraint equation is obtained by varying lapse function A^: 



T,Y'^ + TsX'' + poe 



-2X 



The behavior of the solution depend crucially on the sign of Fi 



Fi < case: 



- cosh cr 
c 



X = In 
Y = At + B. 



where, c, A, B and q = ^/j^ are arbitrary real constants. Using the cosntraint 



tion, we can determine A in terms of other variables 



^^^^/-r3 ^^y3(2^ + 3) 



|ri| 



Fi > case: 



X = In 



sinh cr I 



Y = Ar + B. 



Having been found X and Y, a and can be found from the relation Eq.(^. 

In next section, we will find t(r) and find that its behavior depends on uj and 7. 



III. PHASE SPACE CLASSIFICATION IN TERMS OF T AND r 



A. Fi < case 
From Eqs.(^ and (0), t^r) is found to be 



t-t,= / rfrei(^+^)(^^+^) 



- cosh (cr) 
c 



(l+7)-l 



(14) 



It is easy to see t(r) is a monotonic function. When r goes to ±oo, t can be approximately 
integrated to be 

t ~ ^e^^\ (15) 

where 

± 2ri V T /; -r 



l+7) + l 



(16) 



-4ri 

We have fixed c = 1. We will say that t is supermonotonic function of r when it is monotonic 
and t runs the entire real line when r does. When t is supermonotonic function of r, 
the universe evolves from infinite past to infinite future. Otherwise, the universe has a 
starting(ending) point at a finite cosmic time ti(t/). The running range of t depend on the 
sign of T. 

- oo < t < oo if T_ < < T+, 
-oo <t<tf if r_ < and T+ < 0, 
ti<t<oo if r_ > and T+ > 0, 

ti<t<tf if T+ < < T_. (17) 

The solution for r_ < is found to be 

3 4 

— - < Lu < — -, and LU > uJk, (18) 



where we have defined 



(37 - 5) _ ^ 
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To get Eq.(l^) we used following identity. 



3(2^ + 3) = (^)2_^ri. 

The region II in Fig.l is correspond to this solution. For T_ > 0, we have solution: 

3 4 
uj > --, uj <u^, or u > --. (19) 

The region I and VII is satisfy this solution. 
If r+ < 0, the solution to Eq.(|T|) is 

3 4, 

< UJ < — , and UJ < cj^. (20) 

2 3 

The region I correspond to this solution. For > 0, we have 

3 4 
UJ > — , UJ > uj^, or UJ > — . (21) 
2 3 

The region II and VII satisfy this solution. 



B. Ti > case 



In this case the solution X{t) in Eq.(|l3|) has a singularity at r = 0. So we have to look 
at the behavior of t near r = carefully. From Eqs.(^ and (0), the t(r) can be written as 



t-tn 



dre 



|(l+7)[^^^-+B] 



-I sinh(cr) 

c 



-|^(l+7)-l 



The asymptotic behavior of t in the limit r —>■ ±oo, is given by 

1 



where 

3^3(2^ + 3) ^^ 
T± = (1 + 7)T 

Notice the difference from Eq.(|I3|). 
The condition T_ < gives solution: 



4ri 



;i+7) + i 



(22) 



(23) 



(24) 



3 4 
uj > — -, uj < uj^^, and u > —-. (25) 



There is no region satisfying this solution. For TL > the solution is 



u > --, uj> uJn, or UJ < --. (26) 



Therefore III, IV, V and VI satisfies this solution. 
The solution for T+ < is 



3 , 4 

UJ > — -, UJ > uJk, and uj > —-. [27) 

The region V and VI satisfy this solution. For T+ > the solution is 

3 4 
UJ > --, UJ <uj^, or UJ < --. (28) 

The region III and IV satisfy this solution. 

So far our analysis is parallel to the previous section. However, we have to pay attention 
to the behavior of t(r) near r = 0. In the limit r — ^ 0, 

t^^H-. (29) 
1 — r} 

where rj — + 7) + 1. Notice t{r) is regular at r = 0, if 77 < 1. When 77 > 1, t{T) is 

singular at r = 0. So we consider ^(t) in the region —00 < r < and < r < 00 separately. 
The condition 77 > 1 is equivalent to 

UJ > tura- (30) 

When T goes to zero from the below, 

' - -ITT' P^' 

which means t — > 00 as r — > —0. On the other hand, when r goes to zero from the above, 

tr^-^-^l—, (32) 
77-1 

so that t — > —00 as T — > +0. 

From these analysis, we see that the parameter space of 7 and uj is divided into seven 
regions as depicted in Fig.l. 
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FIG. 1. The phases classified by the relation between t and r. 



We summarize the results that are found in this section. 

• Region I: T+ < 0, TL > 0, Fi < 0; i evolves from initial time U to final time tf for 
r e (— oo, oo). 

• Region II: T_ < 0, r+ > 0, Fi < 0; t evolves from negative infinity to positive infinity 
for r e (— oo, oo). 

• Region III: T+ > 0, T_ > 0, Fi > 0; t evolves from initial time ti to positive infinity 
for r e (— oo, oo). 

• Region IV: T_ > 0, T+ > 0, Fi > 0; Since r = is singular, the region of r is divided 
into two regions, t evolves from initial time ti to positive infinity for r e (— oo,0); t 
evolves from negative infinity to positive infinity for r e (0, oo). 

• Region V: < 0, TL > 0, Fi > 0; i evolves from initial time ti to final time tf for 

T e (— oo, oo). 

• Region VI: 7+ < 0, TL > 0, Fi > 0; Since t is singular at r = 0, we should divide 

into two regions, t evolves from initial time to positive infinity for r G (— oo, 0) and 
negative infinity to final time for r e (0, oo). 
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• Region VII: T_ > 0, > 0, Fi < 0; t evolves from initial time ti to positive infinity 
for T G (— oo, oo). 

IV. THE BEHAVIOR OF THE SCALE FACTOR 

We now study the phases of the cosmology by looking at the scale factor a(r) = 
exp(a(r)). 

A. Ti < case 

In this case q;(t) in scale factor e"*-'^-' is given by 



c^3{2uj + 3) 
a r = — r + B 



-Fi 2ri 

In the limit r ±oo, the scale factor can be written as 



In - cosh(cr) 
c 



(33) 



a(r) ~ e^±", (34) 

where the H± is defined by 



= -^^^r T ^c. (35) 

Eg . (1351) for H_ < gives the solution 

u > ujr2, and a; < 0. (36) 

The region II in Fig. 2 satisfies this solution. For > the solution is 

u < ujr^, or u > 0. (37) 

The region I and VI satisfy this solution. 
If if+ < 0, the solution is given by 

uj < ur2, and u < 0. (38) 
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The region I satisfies this solution. For if+ > the solution is 

ijj > U!r2, or cu > 0. 
The region II and VI satisfy this solution. 

B. Ti > case 

In this case, the Q;(r) in scale factor e"*^"^^ is given by 



c^3(2a; + 3) 

air = — r + B 

^ ^ Ti 2ri 

In the limit r —>■ ±oo, a(T) is given by 

a(r) - e^±^ 

where H± is defined by 



In -| sinh(cT) 



c 



c^3(2u; + 3) 
^ Ti ^ 2ri 



The solution to the condition < is 

CO < cur^, cht^d <^ > 0. 

There is no region satisfying this solution. For H_ > the solution is 

cu > u!r2, or lu < 0. 

The satisfying region is III, IV and V. 
The solution for < case is 

cu > U!r2, O'TT'd cu > 0. 

The region V satisfies this solution. For > the solution is given h 

cu < cur^, or cu <0. 
10 



The region III and IV satisfy this solution. 

Now we consider r ^ hmit. In this hmit a(r) is approximately given by 



a r ~ r 



(47) 



Therefore, if a; > UJT2, goes to infinite as r — > 0. Prom these analyses we get a phase 
diagram Fig. 2. 




-0.5 



-0.5 -1/3 

FIG. 2. The parameter space is classified by the scale factor a(r) 
Summarizingly, we have following cases. 

• Region I: > 0, < 0, Fi < 0; The scale factor a(T) goes to zero size as r — > ±00. 

• Region II: if_ < 0, > 0, Fi < 0; a(T) goes to infinity as r — > ±00. 

• Region III: if_ > 0, if+ > 0, Fi > 0; In this region, the behavior of t is not singular, 
so we need not consider the behavior of a(T) at r = where the scale factor vanishes. 
a(r) goes to zero as r — >• —00 and a(r) goes to infinity as r — >• 00. 

• Region IV: H_ > 0, > 0, Fi > 0; The a(r) goes to infinite size as r goes to zero 
for T e (—00, 0) and zero size as r ^ — 00. a(T) goes to infinity as r goes to zero for 
T e (0, 00) and goes to positive infinity as r — > 00. 
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• Region V: if_ > 0, < ,ri > 0; a(T) goes to infinity as r goes to zero and a(T) 
goes to zero as r goes to negative infinity for r e (— oo,0). a(r) goes to infinity as 
Tgoes to zero and a goes tozero r infinity for r e (0, oo). 

• Region VI: > 0, > 0, Fi < 0; In the limit of r — > — oo, a(T) goes to zero size. 
In the limit of r — > oo, a(T) goes to infinite size. 



V. PHASES OF THE COSMOLOGY 



In previous sections we have studied ^(t) and a(r). Prom all considerations of these 
results, we can classify the parameter space of 7 and lu by the behavior of a{t) into sixteen 
phases. In Fig. 3, we show the phase diagram. In asymptotic region where r — > ±00, we can 
write the scale factor a{t) as: 

a{t)^[T4t-U)]^ 

a{t)^[T+{t-tf)]^. (48) 

From these relations we see that the behavior of the scale factor a{t) depends on the sign of 
T± and the value of H±/T± determins acceleration or deceleration of the scale factor which 
is discussed below. In Fig.4 - Fig. 6, we show the behavior of scale factor by numerical study. 
The — sign in V— indicates the branch for r e (—00, 0). Similarly, V+ means r e (0, 00). 
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FIG. 3. The parameter space is cissified by the behavior of a{t) 
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• Region I, T_ > 0, T+ < 0, if_ > and if+ < 0. The universe evolves from zero size at 
a finite initial time ti to a zero size at a finite final time tf. This region contains the 
matter for inflation 7 = — 1. However, we see in the limit r — > — 00 the cosmic time t 
approaches to ti. 

• Region II, TL < 0, T+ > 0, H_ < and > 0. The universe evolves from infinite 
to infinite size as t runs from negative infinity to positive infinity. This is the region 
where we can find the non-singular behavior of a{t). 

• Region III, T+ > 0, T_ > 0, H_ > and H_^_ > 0. The universe evolves from zero size 
to infinite as t runs from finite initial time to infinity for r e (—00,00). During 
evolution the universe goes to zero as r goes to zero. So we divided into two branches 
at T = 0. 

• Region IV, T_ > 0, T+ > 0, iJ_ > and > 0. The universe evolves from zero to 
infinite size as t runs from finite ti to infinity. Like region /// we divided into two 
branches since the universe becomes zero as r goes to zero. 

By numerical study we depicted the behavior of scale factor a{t). 




FIG. 4. The behavior of the scale factor from phase / to IV 

• Region V, T_ > 0, T+ > 0, i7_ > and > 0. The universe evolves from zero to 
infinite size as t runs finite initial time ti to infinity and the universe evolves from 
infinite to infinite as t runs negative infinity to infinity. 

• Region VI, T_ > 0, T+ < 0, H_ > and > 0. The universe evolves from zero to 
zero size as t runs from initial time ti to infinity and the universe evolves from zero to 
infinite size as t runs negative infinity to finite final time tf. 



• Region VII, r_ > 0, r+ > 0, < and > 0. The universe evolves from infinite 
to infinite size as t runs from finite initial time ti to infinity. 




0.5 1 1.5 -20 -0.8 1.2 3.2 

(a) phase V— {b) phase V+ (c) phase VI— 




-2 0.5 1 1.5 2 2.5 

{d) phase VI+ (e) phase VII 

FIG. 5. The behavior of the scale factor from phase V to VII 
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• Region VIII, TL > 0, T+ < 0, if_ > and if+ > 0. The universe evolves from zero 
to infinite size as t runs from finite initial time ti to infinity and the universe evolves 
from infinite to infinite size as t runs negative infinity to finite final time tf. 

• Region IX, T_ > 0,T+ < O.H^ > and < 0. The universe evolves from zero to 
infinite size as t runs from finite initial time to infity and the universe evolves from 
infinite size to zero as t runs from negative infinity to finite final time tf. 

• Region X, T— > 0, 7+ > 0, H_ > and if+ > 0. The universe evolves from zero to 
infinite size as t runs from finite initial time to infinity. 




(a) phase VIII- (b) phase VIII+ (c) phase IX- 




-3 -1 1 -0.001 0.499 0.999 

(d) phase IX+ (e) phase X 

FIG. 6. The behavior of the scale factor from phase VIII to X 
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VI. ACCELERATION / DECELERATION PHASE 



Notice that as we have seen in Eq.(59) not only the sign of H±/T± but also that of 
H±IT± — 1 is important because the universe will accelerate or deccelerate according to the 
sign of the latter. 



A. Ti < case 

1. H_/T_ > 1 

For T_ > 0, the condition H^/T^ > 1 is reduced to 

(37 + l)v'3(2cu + 3) < 3(1 - 7)(2cu + 3). (49) 

Consider 7 > — | case first. To our surprise, the inequality (^9]) gives us > '^^^'!Zf)l = ^^ri, 
namely Fi > 0. This is contradiction. Now for 7 < — |, Eq.(59) gives u < cupi- The region 
I in figure 3 corresponds to this case. 

For T_ < 0, the H_/T_ > 1 is reduced to 



(37 + l)V3(2cu + 3) > 3(1 - 7)(2cu + 3). (50) 

This inequality gives u < ur-^ for 7 > — |. The region II corresponds to this. For 7 < —1/3, 
we get the condition u > uj^-^ which contradicts to Fi < 0. In a summary, region I and II 
satisfies i/_/T_ > 1, Fi < 0. 



2. H+/T+ > 1 

For T+ > 0, the condition H^/T^ > 1 is reduced to 

(37 + 1)^3{2lu + 3) < -3(1 - 7)(2c^ + 3). (51) 

The analysis is completely similar to the case A. For 7 < —1/3, from above inequality we 
get u < cufi . Since T+ > is satisfied only by the region II, III, IV, VII, it is easy to 
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see that there is no region satisfying all three conditions, Fi < 0,T+ > 0,7 < —1/3. For 
7 > —1/3 Eg. ([51]) gives uj > uj^^ which contradicts to Fi < 0. 
For T+ < 0, H^/T^ > 1 is reduced to 

(37 + l)v/3(2cj + 3) > -3(1 - 7)(2cj + 3). (52) 

When 7 > —1/3, the above inequality gives u < oj^^- Since the condition < is only 
satisfied by the region I, V, VI, there is no region satisfying three conditions Fi < 0, T+ < 
0,7 > —1/3. When 7 < —1/3, Eq. ([5l|) gives 7 > c^n which contracdicts to Fi < 0. In a 
summary, it is always //+/T+ < 1, Fi < 0. 

B. Fi > case 

1. H_/T_ > 1 

For T_ > 0, the condition H_/T_ > 1 is reduced to 

(37 + 1)^3(20; + 3) < -3(1 - 7)(2u; + 3). (53) 

Consider 7 > —1/3 case first. The left hand side is always positive while the right hand 
side is always negative. So there is no solution for this. For 7 < —1/3 case. Eq.(p3D gives 
UJ < cupi which contradicts to Fi > 0. 

For T_ < case, from section III, we know that there is no region satisfying T_ < 
condition. Therefore there are no solutions for conditions Fi > 0,H^/T^ > 1. We 
summarize if Fi > then we have H^/T^ < 1 for all region. 

2. H+/T+ > 1 

For T+ > 0, H^/T^ > 1 is reduced to 

(37 + l)v^3(2cu + 3) < 3(1 - 7)(2cu + 3). (54) 
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Consider 7 > —1/3 case. The above inequality gives uj > c<Jri- Part of the region V satisfies 
these conditions. For the solution to 7 < —1/3 case, Eq.(|^) gives u < which contradicts 
to Ti > 0. 

For T+ < 0, H^/T^ > 1 is reduced to 

(37 + l)^j3{2uj + 3) > 3(1 - 7)(2c^ + 3). (55) 

When 7 < —1/3, we find solution uj > c^ri which is satisfied by the region VI in Fig 3. For 
7 > —1/3 case, Eq . (|55D gives u < ur-^ which contradicts Fi > 0. 
We summarize what we have obtained so far by a table. 



phase 


sign of 
Fi 


sign of 


sign of 


range of 
t 


H^/T_ 




I 


— 


+ 


— 




H^/T^ > 1 


< H+/T+ < 1 


II 


— 


— 


+ 


(—00, 00) 


H_/T_ > 1 


< H+/T+ < 1 


III- 


+ 


+ 






< i/_/r_ < 1 


> 1 


III+ 


+ 




+ 


[ti, 00] 






IV- 


+ 


+ 




[tutf] 


< H_/T_ < 1 


H+/T^ > 1 


IV+ 


+ 




+ 


[ti, 00) 






V- 


+ 


+ 




[ti, 00) 


< iJ_/T_ < 1 


H+/T+ > 1 


V+ 


+ 




+ 


(—00, 00) 






VI- 


+ 


+ 




[ti,tf] 


< H_/T_ < 1 


H+/T^ < 


VI+ 


+ 






[ti, 00) 






VII 




+ 


+ 


[ti, 00) 


H^/T^ < 


< H+/T+ < 1 


VIII- 


+ 


+ 




[ti, 00) 


< H_/T_ < 1 


H+/T^ < 


VIII+ 


+ 






{-00, tf] 






IX- 


+ 


+ 




[ti, 00) 


< iJ_/r_ < 1 


< H+/T+ < 1 


IX+ 


+ 






(-00, t/] 






X 




+ 


+ 


[ti, 00) 


< i/_/r_ < 1 


< H+/T+ < 1 
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VII. EINSTEIN FRAME 



In this section, we study the cosmology in Einstein frame. Especially, we investigate the 
difference between the behavior of the scale factor in Einstein frame and that in the string 
frame as well as the possibity to avoid the initial singularity. 

The metric in Einstein frame is obtained from string frame metric by transformation 
using the relation qe^u = ^''^diiu- 



dSE^ 



-dts'^ + e^^^dxidx' {i = 1, 2, 3). 



(56) 



From above relations, we see a^; = a — | and dtE = e~2dt. Then we can obtain solutions in 
Einstein frame combining the above relation with the solutions in string frame. Therefore 
the solutions in Einstein frame are 



(37 + 1)^3(2^ + 3) 
C(e{t) = c r — m 



-4ri 



cosh(cr) 



for Fi < 0, and 



(37 + l)v/3(2^7+3) 
<^e{t) = c r - In 



-I sinh(cr) 
c 



'iSl±lEl±^y?2±ls, ,57, 



(3,+ l)r. + 4r,\ 32 + i^_ (58) 



for Fi > 0. In next section we find out ^^(t) and see that the interval of t_E(T) can be 
classified by lj and 7. 

A. Classification of the phases by t^; and r 

1. Fi < case 



From Eq.(^ in section II, and using Eqs. (|56[) and (0), ^^(t) can be read 



dre' 



3(37+1)^3(2^+3) 



3[4ri + (37+l)r2] 



[cosh(cr)] 



(59) 
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In the limit r ±00, we can write 

1 

tE ~ — e 
■J-± 

where 



Te±t 



= 4f^ ^ 8f^ • 

The condition for Te- < is reduced to 



(37 + l)V3(2cu + 3) < 3(1 - 7)(2cu + 3). (61) 

Consider first 7 > — | case. The solution for Eg .(1611) is > Ugo which violate Fi < 0. For 
7 < — we get uj < cupi- Now Te- > case that is 



(37 + l)V3(2cu + 3) > 3(1 - 7)(2u; + 3). (62) 

For 7 > — i, Eq.(p^) gives u < c^n- However, for 7 < — |, the left hand side is negative 
while right hand side is positive which is inconsistent. 
Now consider Te+ < case which is reduced to 



(37 + l)v'3(2cj + 3) < -3(1 - 7)(2cu + 3). (63) 

For 7 > — |, in Eq. (|63|) the left hand side is positive while right hand side negative which is 
inconsitent. For 7 < — |, we have uj < cx^ri- Consider Te+ which is reduced to 



(37 + 1)^3(2^ + 3) > -3(1 - 7)(2u; + 3). (64) 

For 7 > — |, we have u < cjri- For 7 < — |, Eq.(|6^) gives u > ujr-^ which violates Fi < 0. 
One can summarize that under the condition uj < cjpn 

For 7 < : Te- < 0, Te+ < 0. 

For 7 > : Te- > 0, Te+ > 0. 
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These relations tell us that 

For 7 < : -oo < ti? < t^/. 

For 7 > : tEi<tE < oo. 

We emphasize that there is no region where Ie can run from — oo to +00. This is because 
Te+ and Te- have the same sign in any given region unlike the string frame. 



Ti > case 

In this case ^^(t) can be read from Eqs.@, (|56|) and (^8]), 
tE -tEo = I c/re3°^(-) 

3(37+l)V3(2i^+3) 3(4ri+(37+l)r2) 

dre sinh(cr)| _ (65) 

As we saw above, this case is singular as r — » 0. So it is necessary to consider the behavior 
in that case. As r — 0, 

tij~sign(r)- , (66) 

1 — rj 

where r] = ^l(^'i+'^)^2+4:ri] ^ ^ ^ 1, t^; is singular at r = 0, while if 77 < 1, it is regular, rj > 1 
case gives 

(5 + 37) 
3(1 -7^) 

The other case, f] < 1, gives < cu* which does not overlap with Fi > 0. Therefore there is 
no regular region for Fi > 0. As r ^ —00, +00 while as r +0, Ie —00. Let us 

find out the behavior in the region r —>■ ±00. In this limit tE and r is given by 



tE — tEO 



J dTe^^^\ (67) 



where 



3(37 + l)V3(2^ + 3) 3[(37 + l)F2 + 4Fi] 

4f^ ^ sf;^ • (^^^ 
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The condition Te- > is reduced to 

(37 + l)^?,{2uo + 3) > -3(1 - 7)(2cu + 3). (69) 

For 7 > — |, the left hand side is always positive while the right hand side always negative. 
So Eq.(^) satisfies trivially. For 7 < — |, Eq.(^) gives ui > ur-^ which is consistent with 
Fi > 0. 

For Te- < 0, we have the inequality 

(37 + 1) V'3(2cj + 3) < -3(1 - 7)(2cu + 3). (70) 

For -y > since after dividing by 37 + 1 the left hand side is always positive while the right 
hand side is always negative. So there is no solution. For 7 < — |, Eq. (|70D gives u < ur-^ 
which contradicts to Fi > 0. 

Now consider the case Te+ > 0. This condition is reduced to 

(37 + l)^j3{2uj + 3) > 3(1 - 7)(2c^ + 3). (71) 

Consider first 7 > — |. Eq.(|7TD gives u < Ur-^ which contradicts to Fi > 0. For 7 < — |, the 
left hand side is always positive while the right hand side always is negative. So we have no 
solution for Te+ > 0. 

Consider Te+ < which is reduced to 

(37 + l)73(2cu + 3) < 3(1 - 7)(2cu + 3). (72) 

For 7 > — |, Eq.([7^) gives u > ur-^ which satisfies Fi > 0. For 7 < — |, in Eq. ([7^) after 
dividing by 37 + 1 we see that the left hand side is positive while the right hand side is 
negative. It is easy to check that when uj > cjp^, 

Te- < 0, and Te+ < 0. 

In a summary, for Fi > 0, t^; is singular as r 0. runs from initial time tEi to +00 for 
r G {—00, 0) and —00 to final time tEf for r G (0, 00). 
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B. The scale factor 

In previous section we see that ^^(t) is not monotonia function which is crucial to decide 
whether the scale factor is singular or not. In this section following the same procedure we 
study the behavior of the scale factor, a_B(r), and classify by uo and 7. 



1. Fi < case 



For this case we use the Eq.(^), then 



aE[T) = e 



, (37+l)V3(2^+3) _ 



cie 



- cosh(cr) 
c 



(37+i)r2+4ri 



(37+l)B 

where Ci = e 4 



In the limit r ±00, the scale factor can be rewritten as e^'^^'^ where 



(37 + 1)^3(20; + 3) (37 + 1)12 + 411 
^E± = 77; T 



(73) 



(74) 



4ri SFi 

Since, in the limit r ±00, He and Te are proportional(TE = 3He), we can analyse by 
using Te in the previous subsection. 
Summary: 

aE(r) runs from + 00 to for 7 < — -. 

O 

a£;(r) runs from to +00 forj > —-. 

o 



2. Ti> case 



From Eq.(ES|), we can write the scale factor for this case as follows. 



aE[T) = e 
= 026 



(37+l)V3(2^+3) [(37+l)r2+4ri] 

*ri |sinh(cr)| 



(75) 
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(3-y+l)B [(37+l)r2+4ri] 

where C2 = e 4 (2j st^ _ Since as we saw in the section VII. l.B, rj < 1 and 
Fi > are not consitent with each other, we consider the case rj > 1. In this region there is 
singularity at r = 0. So we need to consider the limit t —>■ 0: 



-v/3 
(37+i)r2+4ri 



sr, 



(76) 



As r — > 0, the behavior of the scale factor always goes to positive infnity because 77 > 1. As 
r — i> ±00, a£;(r) can be written as 



where 



(37 + 1)^3(20; + 3) [(37 + m + 4ri] 
= 4f^ ^ 8f^ ■ ^^^^ 

By the same analysis of subsection VII. l.B we can write the solutions. Under the condition 

He- > and He+ < 0. 

The scale factor evolves from to +00 for r G {—00, 0) and from 00 to for r G (0, 00). 

Now we study aE{tE)- First consider in the limit r ^ 0. In this limit Fi < case is 
regular. Therefore we investigate Fi > case. From Eqs.(|5BD and ( [75| ) we write aEitE) as 

aE{tE) ~ [(1 - r])sign{T)tE\'^^ ■ (78) 

When 3(;pYj > 1 i-e. 1 < < |, the scale factor will accelerate while when ^ ^ 

?7 > I the scale factor will decelerate. To accelerate, ui should satisfy 

uj > ujr2 for 7 < -3- 
1 

UJ < ujr2 for 7 > -g- 

In Fig. 7, region IV correspond to these relations. For deceleration uj should satisfy 

. 1 
UJ > ujr^ for 7 > --. 
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, 1 

u < ur^ for 7 < --. 

Region III and V satisfy these relations. 

In both cases(ri > or Fi < 0), in the limit r ±oo the scale factor aEitE) behaves 

as 



(79) 



Since we already know that Te and He satisfy Te± = SHe±, from Eq. ([79|) , aEitE) can be 
written 

aEitE) ~ tE^^^. 

It is interesting to compare with Einstein general relativity where aEitE) behaves as 

aEitE) ^ tE^^^ for dust (7 = 0). 
aEitE) ~ tE^^"^ for radiation (7 = 1/3). 



0.5 



-0.5 



-4/3 
-1.5 



IV 




FIG. 7. The classification of phase space in Einstein frame 
We summarize the behavior of the scale factor in terms of t^. 

• Region I: TE-i= SHe-) < 0,TE+i= SHe+) < 0. The universe described by aEitE) 
evolves from infinite size to zero as time runs from negative infinity to finite final time 

tEf- 
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• Region II: Te- > 0,Te+ > 0. The universe evolves from zero size to infinite one as 
time runs from finite initial time tEi to positive infinity. 

• Region III, IV and V: Te- > and Te+ < 0. The universe evolves from zero to infite 
size as time runs from finite initial time tEi to infinity for r G (— oo, 0) and infinte size 
to zero as time runs from negative infinity to finite final time Iej for r G (0, +oo). 
Furthermore, in the limit r ^ the universe has acceleration or deceleration regime 
which depends on the range of rj. Region III and V is decelerationary while region IV 
is accelerationary phase. 



VIII. DISCUSSION AND CONCLUSION 

We have considered string motivated Brans-Dicke(BD) cosmology with perfect fluid type 
matter which arise when a certain kind of the dilaton coupled p-brane gas is dominating the 
universe. This is the complementary study to our earlier work where p-brane gas that 
does not couple to the dilaton was studied. Cosmology is classified into 16 phases according 
to the asymptotic behavior of the time interval and the scale factor. This is qualitatively 
similar to the result obtained before for the dilaton coupled case 0J^. In string frame, there 
is a phase where the cosmology has no singularity, namely, region II in Fig 3. In Einstein 
frame, contrary to the string frame, there is no singularity free phase. This is partly due to 
the difference of cosmic time(t in string frame and t_B in Einstein frame) and partly due to 
the dilaton factor relating two frames. In asymptotic regime, r ±oo, the behavior of the 
scale factor is Ie^^^ ■ To our surprise the inflationary regime of the dilaton-graviton string 
cosmology is gone in the presence of the matter. The matter contribution seems to give 
mass term or potential to the dilaton regulating the dilaton from growing dilaton kinetic 
energy g. 

In a recent study [|T^, with the assumption of holographic principle, it was argued that 
this principle requires the existence of graceful exit by smoothly connecting the pre and post 
big-bang branches. According to the [11|, all cosmological solutions of p-brane dominating 
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the universe can be mapped to the present case or the case studied in In any case, 

there is no solution which exhibit both inflation and graceful exit. Therefore our result 
draws a negative conclusion to what has been said in |iy . 
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